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Abstract. - We deform a two-dimensional (2D) foam, created in a Langmuir monolayer, by 
applying a mechanical perturbation, and simultaneously image it by Brewster angle microscopy. 
We determine the foam stress tensor (through a determination of the 2D gas- liquid line tension, 
2.35 ± 0.4 pj-m -1 ) and the statistical strain tensor, by analyzing the images of the deformed 
structure. We deduce the 2D shear modulus of the foam, fi = 38±3 nN-m -1 . The foam effective 
rigidity is predicted to be 35 ± 3 nN • m _1 , which agrees with the value 37.5 ± 0.8 nN • m _1 
obtained in an independent mechanical measurement. 



Introduction. - A liquid foam, made of polyhedral gas bubbles separated by thin liquid 
walls forming a connected network jj, is a mixture of two fluids. It has nevertheless a solid- 
like elasticity, characterised by a shear modulus /i, proportional to the surface tension of the 
walls [21 E]. In fact, shearing a foam modifies the total length of the walls, thus the foam 
energy. The value of /i can be determined in numerical simulations ^JEHS]; however, it is 
still an open problem to predict analytically its value for a real foam, which has a finite fluid 
fraction and an inherent disorder due to its distribution of bubble sizes. 

Here, we compare two experimental measurements of [i. First, by global mechanical mea- 
surements on the scale of the whole foam, described in terms of elasticity of continuous media. 
Second, and simultaneously, by detailed imaging of the diphasic foam structure, on the local 
level of a few bubbles: this suggests to use two-dimensional (2D) foams. In the literature, 2D 
soap froths have been sheared in Couette geometry, either as bubble rafts El E] or confined 
in Hele-Shaw cells between two parallel plates of glass ^0] ■ 

We investigate the elasticity of a real 2D system: a "Langmuir foam" A monomolec- 
ular layer of amphiphilic molecules deposited at the surface of water ( "Langmuir monolayer" ) 
exhibits a first order transition between a 2D gas phase and a denser 2D liquid (also called 



(*) Author for correspondence at graner@ujf-grcnoble.fr. Fax: (+33) 4 76 63 54 95. 
© EDP Sciences 



2 



EUROPHYSICS LETTERS 



"liquid-expanded") phase. In the 2D gas- liquid coexistence region, the domains spontaneously 
arrange into a foam Walls are stabilised by electrostatic dipolar interactions of the 

molecule themselves, without any added external surfactant |12l ll.ij . Such Langmuir foams 
are approximately characterised by a line tension (see below), hence obey the same Plateau 
rules as other 2D liquid foams |14l IT5] , and display the same generic rheological behavior JHj • 
At the global level, we probe the foam effective rigidity k c fi by measuring the force exerted 
on a rigid obstacle moving relatively to the foam |17| . Independently, at the local level, we 
determine the foam stress tensor (through a determination of the 2D gas-liquid line tension) 
and the statistical strain tensor, by analyzing the bubble deformation (T§|; we deduce the 2D 
shear modulus of the foam |19|. We then compare both measurements, in the frame of linear 
elasticity of continuous isotropic media |2U) . 

Methods. - The 2D Langmuir foam is formed at T = 21°C in a home-made Langmuir 
teflon (PTFE) trough, of dimension a x b — 8 x 2 cm 2 , 0.7 cm deep. Pentadecanoic acid 
(C14H29CO2H) is dissolved in chloroform at a concentration of 3 • 1CP 4 mol-1 . In order 
to prevent the dissociation of acid, the pH of the ultrapure water is set at pH=2 by adding 
hydrochloric acid The air- water interface is cleaned by aspiration, then ~ 10 /A of the 

solution are spread onto the surface of water. After 10 minutes, the chloroform is evaporated, 
and a Langmuir foam of pentadecanoic acid forms. Using teflon compression barriers, we 
choose to adjust the surface fraction of 2D liquid to the limit at which walls are robust and 
clearly visible, which is about 23 %. With an average bubble size ^100 /xm, bubbles are large 
enough for image analysis, and there are enough bubbles to perform statistics. 

The experimental set-up, inspired from ref. |21| . is presented in ref. |17| . We determine 
the resistance exerted by the foam on an obstacle displaced relatively to the foam in the 
horizontal plane. The obstacle is the tip of a vertical rod, actually a denuded glass optical 
fiber (Thorlabs), which also acts as force sensor. Since the shear modulus of the Langmuir 
foam is low, we must here take a much softer fiber than in ref. ^Z]- We chose a length 
L = 3.5 ±0.1 cm, then attack the glass fiber with fluorhydric acid at 40% concentration until 
we reach its core, of diameter 2r = 6 ± 1 /xm. 

The lower end of the fiber plunges vertically in the foam and immerses 10 /jm below the 
surface. The fiber is silanised (n-octadecyltrichlorosilane diluted at 2% in octane): the contact 
angle, measured with a camera attached to the side of the trough, is close to 90°, and the 
residual meniscus is small (white region on Fig. c). 

The upper end of the fiber is held in a concentric chuck fixed on a horizontal translation 
stage coupled to a motor, allowing for a horizontal displacement Xi mp : this applies a horizontal 
point-like deformation to the monolayer. The horizontal deflection £ of the free end of the fiber 
is then measured by connecting the fiber to a laser diode, and collecting (in a photodetector 
placed under the trough) the beam transmitted at the free end of the fiber. We then deduce 
the horizontal resistance force F exerted by the foam on the vertical fiber as F = K(. Here, 
the fiber rigidity K (calibrated by holding the fiber horizontally at one end, and measuring its 
deflection under its own weight) is K = 362 ± 30 nN-m -1 , suitable for precise measurements 
in the pN range. 

The foam effective rigidity k e g, i.e. the ratio of the force F to the horizontal displacement 
X = X imp — C of the fiber within the foam, is measured as fe e g = F/X = K(/X. It 
characterises the whole foam (it depends on its elastic moduli, as well as its size and geometry) 
under a given applied deformation, as long as it remains elastic. It is physically intuitive: k e g 
represents the rigidity one would feel by sticking a finger in the foam and moving it laterally. 

The contrast and resolution of the Brewster angle microscope |23j have been optimised to 
image this foam; for details see ref. |2S] pp. 54-60. Briefly, an incident He-Ne (Uniphase, 30 
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Fig. 1 - Four images of the movie I22|: (a) 10.5 s before, (b) 0.7 s before, (c) 0.7 s after and (d) 
28.4 s after the fiber passes through the middle of the field of view. The fiber moves from the left to 
the right (white arrow) at a constant applied velocity Xim P = 620 /im-s" 1 . The size of each image 
(here corrected by the projection factor at Brewster angle) is 768 x 850 (im 2 . The 2D gas phase 
of the monolayer, of low density, corresponds to a zero reflected intensity and appears black on the 
images, whereas the liquid phase, denser, appears brighter |23|. The average area is larger on the 
last images than on the first ones: this is due to the large-scale inhomogeneities (unavoidable with 
our preparation method) made apparent by the bulk flow, and not to the increase of area of each 
bubble separately (which contributes less than 6 % during our experiment, and is likely due to a 
partial solubilization of amphiphilic molecules in water, rather than to actual coarsening I24|1. We 
do not observe wall breakage (it appears only at 4 times higher velocity, or with a stiffer fiber), nor 
perturbation of the liquid-gas coexistence. 

mW, 632.88 nm) laser beam is tilted at an incident angle equal to the Brewster angle of the 
air- water interface ib ~ 53°, passes a Glan Thomson polariser (Melles-Griot) in the plane of 
incidence, and a quarter wave plate (Melles-Griot). After reflexion on the surface of water, 
the beam enters a combination of two objectives (Zeiss) with small numerical aperture (NA 
= 0.1), magnifying 5 and 10 times, respectively, and an analyser (Newport). The image of 
the surface forms on a CCD camera (sensitivity 10 -3 lux) and is recorded simultaneously on 
video and on computer via a Scion Image frame grabber. 

On each image (Fig. we measure the contribution of the wall network to the stress 
tensor, as follows PS). In a 2D foam, a wall represents two gas-liquid interfaces: thus the wall 
tension is r = 2 A, where A is the gas- liquid line tension (see below). The network contribution 
to the stress a could be measured from the image, by identifying the bubble walls which cross 
a given line of unit length and normal vector v. The vectorial sum of their tensions (20] , here 
t = 2Ae, where e is the unit vector tangent to the wall (its orientation is unimportant in what 
follows), would measure a ■ v. As shown in ref. |18| . we measure the stress tensor with better 
statistics [IS] if we use its equivalent definition over the bulk of the image [27| , which writes 
here (neglecting the curvature of walls [HJ): a = S^ 1 J^r ® £, where S is the area of the 
image, I = £e: I is the length of the wall; the sum is taken over all the walls present on the 
image; <g> is the tensor product: <7y = S^ 1 J2 T i^j CHI- 

We thus need to determine the gas-liquid line tension A (Fig. |2t)- We adapt the method 
of ref. as follows. We record images in fluorescence microscopy, with 1% NBD-HDA 

dye (all other conditions being unchanged). With a localised laser heating (2 K), we break 
one wall. The resorption of the broken wall's free extremity is driven by the wall tension 
r = 2A. It works against the dissipation force, -Fdrag- Since the monolayer is not in a dense 
phase, the usual 3D viscosity 773D of the water subphase dominates the 2D one [321 ESI CS] ) 
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Fig. 2 - Measurement of the line tension A. (a) Image (320 x 215 pro 2 ) by fluorescence microscopy. 
By breaking a wall (on the right of the image, not shown), we obtain a bubble (radius p), attached 
to a single wall (length x), which retracts due to the line tension, (b) Plot of wall length x(t) versus 
J* dt/ p{t) (see text) for four experiments with different initial length x(0). We obtain straight lines, 
which proves that, at these length scales, the line tension is constant despite long-range dipolar 
interactions [T21 1131 l2"31 



and: -Fdrag = — (37r 2 /4) p 773D x. Here x is the rate of retraction of the length x of the 
wall; p is the radius of the object at the end of the wall. We chose the largest possible p (to 
enhance sensitivity): here, a free bubble remains attached at the end of the wall (Fig. The 
prefactor (37r 2 /4) *=s 7.4 is calculated at the limit of vanishing surface viscosity (zero Boussinesq 
number) for a deformable object (here: the bubble) |35| . as opposed to a solid object which 
would lead to a prefactor of 8 [2011^31 ■ We then deduce A from the balance of both forces. Since 
the bubble radius p(t) can vary with time, we write: x(t) — x(0) = — 8A(37t 2 ?73d) _1 Jq dt/p(t). 
From the slopes of Fig. JSJs), we measure A = 2.35 ± 0.4 pj-m -1 . 

To quantify the bubble anisotropy, we use the local texture tensor M |18j . It is a tensor 
constructed using all walls in a local region (here: the field of view) of the foam: M = {t®t), 
i.e. Mij — (li&j)., where I is the vector linking both ends of a wall, and () stands for an average 
over the walls jTHJ ■ Its logarithm InM has the same axes as M, and is real and symmetric 
|36|. We use it to define the statistical strain tensor U = (InM — lnMo)/2, where Mo is the 

reference value of M in the undeformed state of the foam. This tensor U reduces to the usual 
definition of strain in the validity limits of classical elasticity J8| ■ 

We measure the two components U xy and U xx — U yy : they are independent of Mq; in a 
weakly deformed material, they are roughly equal to (i 2 x — l y )/{ll + i^} and (£x£ y )/{£ x + £y), 
respectively. We then compare them to the corresponding components of the stress tensor 
(Fig. |2^) the relation between the stress and strain tensors defines the shear modulus 
p, |2()| . In fact, in a 2D linear, homogeneous and isotropic medium, where the stress has no 
vertical diagonal component ("plane stress", cr zz — 0), the 2D Hooke law defines fj., and the 
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Fig. 3 - Two independent determinations of the Langmuir foam elasticity, (a) Measurement of the 
shear modulus /j, based on analysis of 28 images. Each image provides two points, one for each 
component: closed squares: normal stress difference a xx — a yy versus Uxx — U yy ; open circles: shear 
stress a xy versus U xy . Solid line: linear fit through all points, slope 76 ± 6 nN-m" 1 . (b) Mechanical 
measurement of the effective rigidity k e g- The force F = KC, exerted by the foam on the fiber is plotted 
versus the displacement X = X- lmp — £ of the fiber free end. Closed squares: the fiber is displaced 
in the direction Xj mp > at a constant velocity of 620 firn-s" 1 . Open squares: it returns back to its 
initial position at the same velocity. Solid line: linear fit through all points, slope 37.5 ±0.8 nN-m . 
Open circles: control experiment, at the surface of pure water. 



2D Poisson ratio v (— 1 < v < 1), through |17| : 

® xx ~t~ ®yy ^ ~ {Uxx ~1~ ^yy); ®xx ®yy (J-^xx ^yy)' ^xy 2fiU X y. (1) 

Results. - Initially (Fig. QJl), the bubbles are relaxed and the foam is isotropic. The fiber 
approaches the field of view (Fig. and compresses the bubbles. Behind the fiber (Fig. 
^i), the bubbles are stretched. After relaxation (Fig. ^i) the foam becomes again isotropic. 
Fig. shows a plot of the elastic normal stress difference a xx — a yy versus U xx — U vv , and 
of the shear stress a xy versus U xy . All the data collapse on the same straight line; its slope 
measures 2fi (eq. [TJ: 

^ = 38 ±3 nN-m" 1 . 

This value can be compared to the theoretical computation for a 2D foam with a regular, 
dry honeycomb structure |57| |23 EH of bubble area A: n h = 0.465 x (2X)A-^ 2 . For our 
system, 2A = 4.7 pJ-m -1 , the average bubble area A -1 / 2 = 1.7- 10 4 m -1 , and this theoretical 
expression would give a prediction, /i = 37.2 nN-m" 1 , similar (within errors) to our measure- 
ment from image analysis. A previous experimental measurement on a dry (fluid fraction 
< 3%) 2D soap froth found fi higher than for the honeycomb of 20%, and suggested to inter- 
pret it as an effect of disorder in wall lengths (T^J- Here, we cannot discuss such effect, since 
it would be counterbalanced by the effects of the fluid fraction (our 2D wet foam is expected 
to have a lower (x than the theoretical dry one pQ, similarly to what happens in 3D [51 139p. 

This local measurement of the shear modulus leads to a testable prediction: the value of 
the global foam rigidity, experienced by the free end of the fiber. The foam rigidity k c g is a 
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function of \i and v\ due to the logarithmic range of 2D elasticity |17|. the rigidity also depends 
on the set-up geometry and the boundary conditions at the edges of the trough, which we do 
not know. We use for simplicity the calculations for non-slip boundary conditions [T2|. Here 
the fiber radius is r = 3 /xm, the trough width and length are a — b = 2 cm, and we obtain: 

kcS ~ 1.085 - 0.362 i/' ^ 

The foam Poisson ratio v, which is between —1 and 1, is not measurable in our experiment. 
Since it is probably determined mostly by the network of bubble walls (in the liquid-gas 
coexistence region, the pressure inside the bubbles is constant j!4j , hence does not contribute 
to the foam compressibility), it is likely to be much smaller than unity 0U| : this is compatible 
with the visual impression [25]. If we neglected it for simplicity, eq. J2J would yield, using 
our measured value fi = 38 ± 3 nN • m . 

keg (image) = 35 ± 3 nN • m . 

This value predicted under all the above assumptions can now be compared to the direct 
mechanical measurement of fc e ff, shown on Fig. ^p). The subphase contribution to the 
force is negligible, as shown when performing the same experiment on the surface of pure 
water. During the deformation (Fig. [fj, the foam apparently remains in its elastic regime. In 
fact, first, the force-displacement plot is affine over a displacement of 3 mm, despite the (yet 
unexplained) initial time-lag visible on Fig. (Ph)- Second, as is visible on the movie [221, 110 
bubble rearrangement and no wall breakage is observed. Third, even along the path of the 
fiber, the foam is intact: when the fiber is displaced backwards, the force-displacement plot 
is reversible. We find that its slope k c g agrees (within errors) with the above prediction: 

/feoff (force) = 37.5 ± 0.8 nN • m _1 . 

Perspectives. This result illustrates how adapted is the statistical strain tensor to 
describe locally the elastic properties of disordered media. Beside foams, we intend to adapt 
it to other amorphous materials, such as glasses and polymer networks, to investigate how the 
microscopic disorder of the structure affects the material's elasticity tensor. 
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